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Very precise measurements of masses and cross sections, with errors of < 1%, are expected to be achiev- 
able with a future linear collider. Such an accuracy gives sensitivity at the level of quantum corrections, 
which therefore must be incorporated in order to make meaningful predictions for the underlying new 
physics parameters. For the chargino-neutralino sector, this involves fitting one-loop predictions to 
expected measurements of the cross sections and forward-backward asymmetries and of the accessible 
chargino and neutralino masses. Our analysis shows how an accurate determination of the desired pa- 
rameters is possible, providing in addition access to the mass of the lighter stop. 
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1 Introduction 

A linear collider (LC) will be an ideal environment for high precision studies of physics beyond the standard 
model (BSM). A particularly well-motivated BSM theory is the Minimal Supersymmetric Standard Model 
(MSSM). Fits to the latest LHC data, see e.g. YU, suggest that, if nature can indeed be described in terms 
of the MSSM, charginos and neutralinos could be among the lighter supersymmetric (SUSY) particles, and 
therefore we investigate what the LC can reveal about the structure of the chargino-neutralino sector of the 
MSSM. 

At leading order (LO), the chargino-neutralino sector can be parameterised via the gaugino masses Mi 
and M 2 , the higgsino mass \x and tan (3, the ratio of the vacuum expectation values of the two neutral Higgs 
doublet helds, see e.g. ref. [2] for details. If elements of the chargino and neutralino spectra are indeed within 
reach of a linear collider, the determination of Mi , M 2 , /i and tan (3 at the percent level via pair-production 
has been shown to be possible at LO (see e.g. ref. [3]). However, one- loop effects in the MSSM can be large 
and therefore higher order calculations of processes in the MSSM are crucial for accurate predictions within 
this model. On taking these corrections into account, additional MSSM parameters become relevant, such 
as the masses of the stops which are so far only weakly constrained by the LHC. By fitting experimental 
results to loop corrected predictions, calculated in the on-shell scheme, it should be possible to extract the 
parameters of the chargino-neutralino sector of the MSSM Lagrangian, as well as to indirectly gain insight 
into the parameters of other sectors contributing via loops. 

In sec. [2] we will outline our strategy and introduce the necessary notation. We will then briefly describe 
the calculation of the loop corrections in sec. [3] including details of the renormalisation scheme used. In 
sec. [4] we will further discuss the method employed in order to fit the MSSM parameters, and present our 
results. Finally in sec. [5] we will discuss the implications of these results. 
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2 Strategy and motivation 



The mass matrices for the charginos and neutralinos, X and Y respectively, are defined by 



X = 



M 2 



y/2M w s p 



Y = 



Mi 


-Mzcpsw 
M z sps w 





M 2 
Mzcpcw 
-M z spc w 



-M z cps w 
Mzcpcw 
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/ 



(1) 



Here M 1; M 2 and /i are the bino, wino and higgsino mass parameters respectively; syy/cw are the sin /cos 
of the electro- weak mixing angle 9w] Mw and Mz are the masses of the W and Z bosons, sp/cp are the 
sin / cos of p. The chargino mass matrix X can be diagonalised by the matrices U and V via the bi- unitary 
transform M x ~ + = U* XV\ and the neutralino mass matrix Y, due to its Majorana character, by a single 

unitary matrix N via M % - = N*YN*. 

At tree level, the fit to the polarised cross-sections of only the light charginos and neutralinos, supple- 
mented by the measurement of the masses of the light particles in the spectrum, namely chi 1 , chi 1 and chi 2 , 
was shown in ref . [3] to be sufficient to obtain Mi , M 2 , \i and tan /3 to an accuracy at the percent level for 
the scenario SPSla E]. Studies have also been carried out with unpolarised beams, using, for example, the 
forward-backward asymmetry instead [5jj7] . The incorporation of loop corrections is strongly motivated by 
the well known observation that these can be relatively large in SUSY compared to the expected experimen- 
tal accuracy. Therefore calculating higher orders is required to ensure that the theoretical precision meets 
the high experimental accuracy achievable at the linear collider. Furthermore, additional parameters enter 
the expressions for the loop corrections, and these can be included in the fit. 

In order to assess the potential of the lin- 
ear collider to access the fundamental parame- 
ters of the MSSM when incorporating loop cor- 
rected theory predictions, we need to choose 
a specific scenario of the MSSM. We ensure 
that the proposed scenario, given explicitly 
in tab. [T] satisfies all existing experimental 
constraints, of which the most stringent come 
from: current LHC limits; the measured Dark 
Matter relic density, which we calculate us- 
ing using mi cr Omegas; the branching ratio for 
b — > S7 and the anomalous magnetic moment 
of the muon. Due to the current status of 
direct LHC searches [8|[9], we further require 
heavy first and second generation squarks and 
a heavy gluino, but take into account the fact 
that the bounds on charginos, neutralinos and 

third generation squarks are much weaker. We consider the sleptons and the heavy Higgs doublet to be at 
the TeV scale, such that they have negligible effect on the size of the loop corrections. Indirect limits lead 
us to choosing mixed gaugino higgsino scenarios, favoured by the relic density measurements, and relatively 
high charged Higgs masses, in light of flavour physics constraints. Assuming this scenario, we carry out the 
following: 



Parameter 


Value 


Parameter 


Value 


Mi 


125 GeV 


M 2 


250 GeV 


/" 


180 GeV 


M H + 


1000 GeV 


M 3 


700 GeV 


tan/3 


10 


M ei . 2 


1500 GeV 


M e3 


1500 GeV 


M h 


1500 GeV 


M qi , 2 


1500 GeV 




400 GeV 


Af 


650 GeV 



Table 1: Table of parameters, where M H + is the mass of the 
charged Higgs, M3 is the gluino mass parameter, M(;/ 9 ). and 
M( e / U v represent the left and right handed mass parameters 
for of a slepton/squark of generation i respectively and Af 
is the trilinear coupling for a sfermion /. 



Calculate one-loop corrections to the amplitude for e + e 



xfxi 



C/0 



Calculate one-loop corrections to the masses for the charginos and neutralinos Xj an d X 

By assessing the sensitivity of these corrections to the MSSM parameters, determine which parameters 
are obtainable from the fit. 

Define an appropriate set of observables which can be used to fit these parameters at loop level, and 
calculate MS H function of the observables 
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Figure 1: Tree- level diagrams for the production of charginos xt an d Xi a t the LC. 



• Fit results as a function of the parameters to the observables and determine the uncertainty to which 
they are attainable. 

3 NLO Calculation 

At leading order, neglecting the electron-Higgs coupling, the process a(e + e~ — > XiXi) is described by the 
three diagrams in fig. [T] The transition matrix element can be written, closely following ref. 10 , 



M aP (e+e -> xtXj ) = ~Q a p [«(e + )7 M w a u( e )] [u( Xj )^w p v{xt)\ , (2) 



in terms of the helicity amplitudes Q a( 3, where e is the electric charge, ujl/r — 1/2(1 T 7s)j o. refers to the 
chirality of the e + e _ current, /3 to that of the xtxj current, and summation over a and (3 is implied. The 
helicity amplitudes can further be defined via 

Q RL =C^ + D Z G R C^ Z Q R n=C« H - i + D z Gn(c« Hjz )\ (3) 
for which the required MSSM couplings for the xtxjli xtxJ Z and eu e xt vertices are given by 

XjXj 7 J X t Xj z c w s w \ J n 2 3Z 

C R + ^=-C L + „-JU^V*), C L + -=— , C* +.-=-— V a , (4) 

xtx } z xtxjZ^ h „ c e+ Xt S W y/2 C pM W V ° e+ *i s w 

and Gl, Gr, D z and Dp are defined via 

I - « 2 

r, 2 b W n _ 8 w 

<^L — 1 <~tR — , 

SwC\V C W 

Dz =7^W z > D *=T^4- (5) 

In the above, D z and Dp refer to the propagators of the Z boson and sneutrino (of mass mp), in terms 
of the Mandelstam variables s and t. We can neglect the non-zero Z width for the considered energies. 
The tree-level cross section in the unpolarised case is then obtained by summing over the squared matrix 
elements, 



K 1 / 2 (s,m-+,m--) 

^trcc _ ' 



64tt 2 s 2 



— f dn^2\M a p\ 2 , where k(x, y, z) = (x - y - zf - Ayz. (6) 

a-fi 



Therefore at leading order we see that the SUSY parameters entering the cross section are the masses of 
the charginos, the elements of U and V, and the sneutrino mass. 
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Figure 2: Example one- loop self-energy, vertex and box diagrams (upper row) and counterterms (lower row) 
for the production of charginos xt an d Xi at the LC. 



The one-loop corrections involve self-energy, vertex and box diagrams, examples of which are shown in 
fig. [2] The diagrams are generated and the amplitudes calculated using FeynArts 
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found in ref. |2j. FormCalc 16-18 was then used to calculate the matrix elements and LoopTools 16 



details can be 
to 



perform the necessary loop integrals. The loop integrals are regularised via dimensional reduction [19f|21| 



16 22 . As seen from 



which ensures that SUSY is preserved, via the implementation described in Refs. 
the diagrams, squarks, sleptons, Higgs particles enter the loops, such that the results now depend on many 
MSSM parameters beyond the small subset at tree level. 

Finite results at one-loop order are obtained by adding the counterterm diagrams shown in fig(2j Although 
FeynArts generates these diagrams, expressions for the counterterms which renormalise the couplings defined 
at tree- level in eq. ^ are required as input, and therefore we provide expressions in terms of the relevant 
renormalisation constants (RCs) for these here explicitly. For the jxfxj > ZxtxJ and ev K x7 vertices, these 
can be expressed as follows, 



XjXj 7 



*i z c w 



sz e . 
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(7) 



where the analogous right-handed parts are obtained by the replacement L — > R, and 
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Note that we define the renormalisation of the chargino fields in the most general way, making use of separate 

L I Ft ~ L I R 

RCs for the incoming and outgoing fields, i.e. coefficients 5Z, / ii and Z ± ' i - respectively for left and right- 
handed charginos, where definitions can be found in ref. [23] . Note that although separate RCs are not 
necessary in the CP conserving case, see the discussion in ref. 24 , we follow this approach to allow the 



analysis to be easily extendable to include complex parameters. The RCs for the SM fields and parameters 
are also defined in ref. 



24 . In addition, we take the counterterm for the sneutrino self energy to be 



5C r ,.r,. 



(SZ^+SZt) P 2 -8ml 



^(SZ Si + 6Zi) 



(8) 



Here SZ„. is the sneutrino's field RC and Sm^. its mass RC, again defined in ref. 23 



for v l = z> e , 

As we renormalise on-shell, the counter terms are obtained by requiring the standard on-shell conditions 
for the mass and field RCs. In order to calculate the mass corrections of the charginos and neutralinos, we 
must renormalise the mass- mixing matrices X and Y via X — > X + 5X and Y — > Y + 5Y. Here SX and SY 
are defined in terms of S Mi, 8M2 and 8/j,. Expressions for these three counter-terms are obtained, following 
e.g. ref. 23 , by choosing three out of the total six physical masses of the charginos and neutralinos to be 
on-shell. For consistency with the Higgs sector, tan /3 is renormalised in the DR scheme. As the process 
we consider has external charginos, and we prefer these to be on-shell, we adopt the NCC scheme which is 
defined in ref. 
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such that Xii xt an d X2 are on-shell 
neutralinos will be shifted, i.e. obtain loop corrections. The expressions for the chargino field RCs can be 
found in ref. 23-25 , where one can also find expressions for the on-shell field and mass RCs for the sneutrino, 



23 



20 



Note that therefore the heavier three 



electron, and gauge bosons, as well as the procedure followed for electric charge renormalisation. Using this 
prescription to renormalise the vertices we obtain UV-finite results. 

Soft and collinear radiation must be included to obtain a result free of infra-red and collinear sin- 
gularities. In the regions E < AE and 9 < AO where AE and AO denote the cut-offs, the radiative 
cross-section can be factorised into analytically integrable expressions proportional to the tree-level cross- 
section er trcc (e + e~ — > xtxj)- The soft contribution can easily be incorporated using FormCalc, however 
the collinear contribution must be added explicitly. At leading order the result for the collinear contribution 
takes the form 



27 28 



J da con ( Pu p2,a-,a + ) = -^ °dx^/ Q (A0)( J da t ^(x Pll p 2l aa^,a + ) + J da tIee ( Pl ,xp 2 ,a-,aa + )) 



where we define the structure functions f + (AQ) 



1 



log 



sAO* 



, /_(A6>) = 1 — x and the limit of 



1 — x y Am 2 

integration xo in terms of AE, xo — 1 — lAEj^fs. Here pi, p 2 and c_, ct + are the momenta and helicities 
of the electron and positron respectively. However, on adding these contributions, the result is cut-off 
dependent (i.e. on AE and AO), and removing this dependence requires a calculation of the full cross section 
for the three body final state, excluding the soft and collinear regions, which we perform using FeynArts 
and FormCalc. 



4 Obtaining fundamental MSSM parameters from the fit 

In the chosen scenario (tab. [T]), the quantities are found to be sensitive to the following: the parameters 
describing the chargino and neutralino sector, namely Mi, M2, fi and tan/3; the electron sneutrino mass 
mp e , as it enters the process at tree level; the stop masses and mixing angle, as the loop corrections depend 
strongly on the stop sector. As observables we include the polarised cross-sections, the forward-backward 
asymmetry (both at beam energies of 350 and 500 GeV) and the masses of the charginos and first three 
neutralinos. 

We summarise the results for the central values of these input observables in our chosen scenario in tab. [2] 



The errors have been estimated using results from Ref. 29 , these are statistical errors only, and no theoretical 



errors are taken into account. We consider the errors both on obtaining masses from the continuum and 
by a more precise threshold scan. The chosen observables are calculated at 1-loop as a function of each of 
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Observable 


Energy/ GeV 


Tree-level value 


Loop correction 


Error 


m~±/GeV 




149.6 


OS 


0.1(0.2) 


m-±/GeV 




292.3 


OS 


0.5(1.0) 


m^o/GeV 


— 


106.9 


OS 


0.2(0.2) 


m^a 1 GeV 




164.0 


2.0 


0.5(1.0) 


m^o/GeV 




188.6 


-1.5 


0.5(1.0) 


o{xiXi) + l& 


350 


2347.5 


-291.3 


1.3 




350 


224.4 


7.6 


0.4 


A FB 


350 


-2.2% 


6.8% 


0.8% 


a(^xr)+/fb 


500 


1450.6 


-24.4 


1 




500 


154.8 


12.7 


0.3 


A FB 


500 


-2.6% 


5.3% 


1% 



Table 2: Input observables for the fit. The errors on the cross sections and forward-backward asymmetries 
are purely statistical obtained for an integrated luminosity C — 200 fb _1 and selection efficiency of e = 0.15. 
The errors on the masses are obtained using threshold scans (continuum). The subscripts + and — indicate 
the beam polarisation configurations (—0.8,0.6) and (0.8, —0.6), respectively, and OS indicates on-shell. 



Parameter 


Loop-level fit (continuum) 


Loop-level fit (threshold scan) 


Tree- level fit 


Mi / GeV 
M 2 /GeV 
/x/GeV 

tan/? 
mp/GcV 

cos &i 
m fl /GeV 
m t ~ 2 /GcV 


125 ±0.4 (±0.8) 
250 ± 1.2 (±2.4) 
180 ±0.4 (±0.8) 
10 ±0.8 (±\\) 
1500 ± 18 (±11) 

sm +220 /+540\ 
8UU_ 170 U280J 


125 ±0.3 (±0.7) 
250 ±0.6 (±1.3) 
180 ±0.4 (±0.8) 
10 ±0.5 (±1) 
1500 ± 24 (± < $ ) ) 

0±0.15 (±8;|) 

fat limits 
wu -120 Ut limit ) 

snn+ 300 /-+ 1000s ! 
»uu_ 170 U290 ) 


122.0 ±0.5 
260.7 ± 1.4 
176.5 ±0.5 

27 ±9 
2230 ± 50 



Table 3: Fit results for masses measured in the continuum and the threshold scan. The numbers in brackets 
are 2a statistical errors for an integrated luminosity C = 200 fb~ . 'Tree level' indicated values calculated 
at tree level are fitted to one- loop corrected observables. Note that the central values of cos 8^, and 
are 0, 400 GeV and 800 GeV respectively. See text for details. 
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the fit parameters, Mi, M%, 
Minuit 



30 



/x, tan/3, m£; e , mj , mj , and We perform a multi-dimensional x 2 fit using 



E 



SOi 



with 



Oi = O i (Mi,M 2 ,fJ,,taxiP,mi> e ,rn il ,m i3 ,0{), 



(9) 



where the sum runs over the input observables Oi with their corresponding experimental uncertainties SOi. 
The theoretical values Oi are calculated at NLO as a function of fitted parameters. The results of the 
fits are shown in tab. [3j We compare the fit accuracy for measurements of the masses obtained from the 
continuum and using threshold scans. Clearly, the more accurate mass measurement allows for a more precise 
determination of the underlying SUSY parameters. But more importantly, it even enables determination 
of further parameters in the stop sector: and 6$, which are otherwise not accessible. Note that here 
the error on mi increases becasue for the fit with masses obtained from the continuum these parameters 
were fixed to their central values. We also show that if values calculated at tree level are fitted to one-loop 
corrected observables, the parameters are not correctly extracted, see by the column labelled 'tree-level'. 
This stresses the importance of including higher order corrections in the analysis. 



5 Conclusion 

We have presented a new NLO analysis showing the possibility to determine the fundamental MSSM pa- 
rameters in the chargino and neutralino sector via chargino production at a linear collider. We calculate 
one-loop corrections to the cross-section, forward backward asymmetry and corrections to the chargino and 
neutralino masses, renormalising in the on-shell scheme for the scenario given in tab.[l] We have summarised 
the central values with statistical errors for an integrated luminosity C = 200 fb _1 in tab. [2} We find that the 
parameters M±, M2, /x, tan/3, mj 2 , mj 2 and cos#t can be extracted from the fit, to good accuracy, as shown 
in tab. [3j We additionally show the crucial role played by an improved measurement of the chargino and 
neutralino masses via threshold scans. Note that this analysis has been extended to include the observables 
rrih and BR(6 — > 57) for a scenario compatible with m/,= 125 GcV in ref. [2]. In addition, we would like to 
incorporate the results of neutralino production, and also investigate the sensitivity to complex parameters 
e.g. the phase of the trilinear coupling of top/stops, 4>A t - 
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